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Abstract 


An underetandlng of the behavior of materials, of dielectrics in particular, 
under Charged particle bombardment is essential to the prediction and prevention 
of the adverse effects of Spacecraft charging. This paper presents an effort to 
Obtain Such an understanding through a combined analytical and experimental 
approach. 

A one-dimensional model for Charging of samples in the LeRc test facility is 
used in conjunction with experimental data taken in this facility to develop "material 
chargttig characteristics" for silvered Teflon. These characteristics are then used 
m a one dimensional model for charging in space to examine expected response. 
Relative charging rates as well as relative charging levels fOr silvered Teflon and 
metal are discussed. 


1. INTRODtJCTIO>i 


Two previous papers*'^ described the test facility, test methods and measure- 
ments, and the results of Various materials characteristics test performed at tiie 
Lewis Research Center in support of the Spacecraft charging investigation. The 
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prortnt papef .umwafUei the ailaljtlcal work whUll has bean pertbrmed Inter- 
actlv(»ly M^lth this exp^flrrierttal work, the goal of the analysis Is t^/vofola. First. 

Is tb model the charging of material samples lii terms of the material's para- 
meters. Second, since a goal of the entire study is to predict behavior of space- 
craft surfaces, an attempt Is made to "Scale" the environment, that is. to relate 
results obtained using a monoenergetlc beam in the ground test facility to expected 
results with distributed particle fluxes of the space environment. An approach to 
establish this environment scaling Is to develop models of charging for both 
charged particle environments, and assume that the material properties are con- 
stant. then differences between material charging behavior under ground test and 
in space are a result of the differences in the two environments. It IS recognized 
that the vacuum levels In the ground test and space environments are also different. 
No attempt is made here to account for. this factor. 

This paper, then, represents a first attempt at attaining me twO goals of 
characterizing material charging and scaling to the space environment. The 
models used ar.i one-dlmenSional and describe charging of samples in terms of the 
charging of a capacitor. This type of model has been used by a number of 
workers^' ^ to describe Spacecraft charging. 

The procedure used herein was to first develop a one- dimensional model to 
describe charging of Sampled in the LeRC vacuum teSt facility. This model con- 
tained a number of oarameters which were varied to provide best fits to experi- 
mental data Obtained in the facility. The values of these parameters which yielded 
the best fit were identified aS the "material charging characteristics. " These were 
then used in conjunction with a ofle-dimefisional model fOr charging Irt the space 
substorm environment to make some predictions of the charging behavior of the^ 
materials in space. The insulator studied here is 5 mil silvered FEP Teflon. 


2. ONK-UIMKNSION \L (l-li) TiKOI’M) TKST MUDKU 

in the LeRC test facility, a monoenergetlc beam of electrons with energies of 
•V. 2 to 20 keV IS directed at normal Incidence to planar samples. A beam current 
density of 1 flA/cm^ was used to obtalri all test data disCuSsed in this report. 

The Ground Test Facility Model Is a quaSlstatlc current balance model, the 
current densities considered are thOSe due to primary (beam) electrons, secondary 
electrons, backscattered electrons, and leakage current through the sample bulk. 
These 'kre denoted by J^. j^. j^g' ^l* respectively. The sample Is assumed to 
charge like a capacitor, thus a time balance equation Is of the form 

dV„ 


4fid 


( 1 ) 



where id the magnitude of the surface voltage. All signs arc explicit in this 
equation, and in all others used in the test facility model, tiiat is. all symbols 
stand for positive numbers. The act-iai surface voltage Is of course negative, and 
this is assumed throughout. Thus this model does not predict positive surface 
voltages correctly, since positive surface voltage would tend to reattract the 
emitted secondary electrons and this effect has not been included here. 

The first term in Eq. (i) represents the net charge deposited on the surface 
<per cm ) in a tin s step, and so is termed the charging current density and denoted 
fay jg. The procedure used was to initialize at t = 0 with Vg = 0. Then was 
calculated from. ® 
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( 2 ) 


The current densities are all functions of and their functional forms are given 

in Figure 1. Equation (2) Is Solved by an iterative procedure. This equation can 
be expressed as 


With th« inltUl MMimpUhn V. . t . 0. a aultable at la chosen, and iV. caloolatad. 

g is then inciemented by AVg (set = AVg for the first iteration) and the procedure 
repeated until equilibrium is reached, that is, until 

dV 

jc= = 0 


The several current densities in the preceding equations are functions of sur- 
face voltage. All but the leakage current density are functions of the primary 
electron beam voltage and current density. The parameters which can be varied 
are the secondary emission maximum yield. and energy for maximum yield, 
V„. the backscatter coefficient. the resistivity, p. and the capacitance C. 

In practice, valued for and were taken from the literature. Values for p 
were determined from the measured surface voltages and currents at equilibrium, 
and fits obtained by varying fand C. Tills is discussed more fully in Section 4. 

It should be noted that this model does not account for beam spread In angle 
or energy, the presence of the vacuum tank walls, or residual gas in the chamber. 
Since it is one-dlmenslonal, it can not, of course, explicitly describe edge effects, 
dr effects due to surface variations. 

Oerlvatiofls of the current density equations are presented in Appendix A. A 
summary of the model and equations is given in Figure 1. 
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Figure 1. 1-D Model Ciround Test Environment 


3. COMPARISON RlTII KXPKRIMKM AL DATA 
3.1 "PToatihg” Aluminum Piute 

The first test data listed in conjunction with the test facility charging model 

were those taken with a bare aluminum platte which is normally used for substrates. 

This plate Was mounted in front of a second identical plate which was grounded to 

the chamber walls. These two plates were held apart and electrically isolated by 

a Teflon spacer plug 0.7 cm long. The plates ate rectangular with dimensions 

15 cm X 20 cm ("“300 cm^ area). Surface voltage of the floating plate was mea» 

sured as a function of time using the surface voltage probe. ^ 

Figure 2 shews these data and the beet fit calculated curves. In obtaining 

these fits to the data, values of 6 _ =3 and = 40o were used. These are 

m max „ 

consistent With ranges of values for these parameters given by Gibbons' for a 
surface layer of AI 2 O 3 . It is reasonable to expect some oxide on the surface of 
the "bare aluminum" plate since it had In fact been exposed to air. This points up 
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Figure 2. Comparison With Experi- 
mental Data 


the strong dept.idence of charging phenomena oh Sample surface condition, and 
indicates that care mrst be taken in making predictions fob charging to consider 
the state of the surface. 

With these values in. the expression for secondary emission, the values of C 
used to obtain the curves in Figure 2 Varied from 1,2 X 10‘^^arads at V 3 = 5 kV 
to 5 X 10 farads at Vg = 16 kV, decreasing approximately linearly with in- 
creasing beam voltage. The values used for the "backScattCr coefficient" varied 
from 0 . 15 at Vg = 5 kV to 0. 5 at Vg = 16 kV, again ih approximately linear 
fashion. The expression used to calculate Secondary emission current density is 
derived from an expression for yield as a function Of primary energy due to 
Sternglass. ® It is plausible that the required variation of 5 to obtain fits to the 
data is accounted for, in part, by deviations of the actual secondary yield from 
that predicted by Sternglass formula. That is, the adjustments to ? represent 
adjustments to the sum of the backscattered and secondary electrons. 

Tne final point of concern here is the time scale for charging of the floating 
pl.ate: it reaches its equilibrium floating voltage With a time constant of Several 
seconds. This iS not surprising, since the capacitance of the test plate to its Sur- 
roundings is expected to be small, llie timescale, is relevant, however, to the 
question of the behavior of composite samples. This is discussed more fully in 
Section 3. 3. 


3.2 SUv«‘r«*d 'IVflon SlunjiU's 

13ie 1-D model Was ne^t used to fit eurrent add SuMace VoU&ge data ffom 
slivered Teflon Sam]ples. These cohsist of an aluminum substrate wim three strips 
of 5 mil silvered Teflon mounted with conduetlve adhesive to the substrate. Each 
of these Strips was 5 cm Wide and 20 cm long. During test, the aluminum substrate 
land consequently the silver) was grounded, while the Teflon Surface was bom*- 
barded with electrons. 

The data and Calculated fits for beam voltages at which equilibrium is reached- 
are shown in Figure 3. These data are a composite of four separate data sets, 
and indicate charging times on the. ^rder of minutes. The error bars reflect the 
Scatter in the data as well as the ±5 percent resolution uncertainty in the voltage 
measurements. Since for insulators there are strong voltage gradients near the 

edges of the samples, the surface Voltage measurements are those read at the 

Sample center which is uniform. 

To obtain these curves, the effective resistance was calculated from the 
equilibrium values of surface voltage and leakage current. These values indicate 
an effective resistivity for the samples Of about 9 X lO^^D-cm, about an order of 
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magnitude below published values tot " eflon. If ohe assumes that the published 
value of 10 tt-em-ls adcurate. theh a parallel path having a resistance ot about 
4 V 10 17 is maleated by the data. This ceuld be a surface leakage, or an edge 

leakage, 6r leakage thrdugh the sheath. 

Values bf 8^ = 3 and = 300 were used in the expression for Secondary 

emission, in aCcotdance with the data given by Willls and Skinner. ^ The values of 
C retjulred ranged frdm 14. 6 pf/cm^ at Vg = 5 kV to 10 pf/cm2 at Vq = 12 kV and 
decreased monotonlcally with- increasing beam voltage. Extrapolation of the curve 
to Vg = 0 indicates a dielectric constant for Teflon of -^2. 1. This decrease in 
effective capacitance is believed to be as sociated with edge effects. Edge gradients 
are observed in the data. T^ese become more pronounced at higher voltages 
reducing the effective area. 2 The value used for the backscatter coefficient 
varies about an order of magnitude for the Teflon samples. Not only does it change 
With beam voltage, but also appears to change during charging at a single beam 
voltage. For the initial portioi of the charging curves. ? varied froto 0. 25 at the 
lower beam voltages to 0. 05 at high beam voltages. At equilibrium, r Varied from 
0. 02 at low beam voltages to 0. 25 at high ones. These variations were not linear; 
rather, t was relatively constant at low and high beam Voltages, with a transition 
occurring between Vg - 8 kV and Vb = 12 kV. Again, part of these variations may 
e due to deviations Of the actual secondary emission from-that calculated by the 
analytical expression being used. 


Some investigation was undertaken to Study the behavior of the Teflon samples 
during arcing. Figure 4(a) Shows a- curve fit to a Composite of two data sets for 
the initial charging transient with a beam voltage of Ifi kV. Figure 4(b) shows the 
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same calculated curved, this ttn;e with a single set of data, th® left hand set of 
cufvee and points are the initial charging tfanslefit. An ahc occurred on this 
sample between the time of- the voltage reading at t = 4 min and that of the current 
reading taken at t = 5-1/2 min; the surface discharged, and Chalglng Was repeated. 
The curves shown for the poet- arc charging transient are Identical to those for the 
initial transient, but Shifted in time. This indicates that the charging transient iS 
quite repeatable, not only from test to test of Teflon Samples, but also through 
arcing. That is, at least short term, the arc ing does not affect the charging 
behavior of Teflon samples. 


:l.:i Suhiplr 

One of the stated objectives of the present work was to investigate the ' environ- 
ment scaling" effects between the ground test and space environments. An obvious 
difference other than the environment between ground tests so far described and 
th’* space condition is that the Studies of slivered Teflon diScusSed above were aU 
conducted with the substrate grounded. In contrast, for the case of a spacecraft 
in orbit, the entire body, including the "grounds" must come into equilibrium with 
the charged particle environment. The question of the behavior of a composite 
system becomes particularly interesting because of tiie divergent time Scales in 
which charging of "noatlrtg" metal plateS and silvered Teflon abOve ground are 
dbsei^ed tb Occur. 

To investigate this question, a composite sample was built and tested. The 
sample consisted of an. electrically floating standard aluminum substrate with two 
Strips of the 5 cm silvered tape mounted on it. Tliis allowed for a 5 cm strip of 
the aluminum between the two Teflon strips to be ejtposed to the electron beam. 

The aluminum substrate was mounted in the test chamber in the same manner aS 
the floating aluminum plate described in Section 3. 1. This configuration iS 
depicted at the top of Figure 5(a). The bottom of this figure shO i r. a voltage trace 
at equilibrium for this sample. 

The eiipectation WaS that this composite Sample would charge In two stages 
because of the different effective capacitances through which the aluminum and the 
Teflon must charge. The prediction, shown in Figure 5(b), IS ba.?ed on the idea 
that when the beam is turned on, the aluminum should charge to Its equilibrium 
voltage with Its time constant of seconds, carrying the teflon voltage with it. 

When this has occurred, the Teflon should continue to charge from the equilibrium 
voltage of the aluminum to its own equilibrium Voltage with its own time constant, 
that is, minutes. The curves in Figure 5(b) were thus obtained by superimposing 
the curves for aluminum above and for Tefloh alonC with the Teflon curve shifted 
so that it coincides With the aluminum curve at the point of equilibration for the 
aluminum. As can be seen from the data plotted in Figure 5(b), the expected be- 
havidr was found. 
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The expectation thotthe aluitiinum and Ttfloh comprising the composite 
sample Should charge to the Same surface voltages as h.d the floating plate: and 
the Teflon samples above grounded substrates Was based On the observation that in 
both those cases the equilibration WaS dominated by secondary emission phenomena 
rather than by leakage currents to ground. TlUs is evidenced in tw« ways. First, 
plots of surface voltage at equilibrium versus beam voltage are straight lines; 
such behavior iS Supposed to be associated with emission dominated equilibration. 
Second, examination of printouts of the model calculations reveals that, at equi- 
librium, the leakage current density term is several Orders of magnitude smaller 
than the othei currents in the model. The conclusion, then, is that for this type 
of composite sample, each part responds to the charging environment with its 
characteristic time constant, and comes into equilibrium at its characteristic 
Surface voltage so long as leakage current does not play a dominant role in the 
equiltbrai ion. 

ft should be noted here that the tests run On this composite sample were not 
extensive. Further experimental investigation of this and other composite samples 
are planned. 
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In order to predict charging behavior In space, a one -dimensional model for 
charging In Space analogous to the one-dimenslonal ground test model waS deve- 
loped. The essential difference between the two models Is that the space model 
assvmes an IsotrOplc MaxwelUan particle distribution containing both electrons and 
ions and a spherical collection geomCtt-y. The current densities are derived from 
Langmuir probe calculations. This type of calculation has been used by several 
authors^' to treat the spacecraft charging problem. Derivation of the 

current density equations is given in Appendix B; a summary description IS pre- 
sented in Figure 6. As Is Indicated ih this figure, the model as presented and 
used here assumes a geomagnetic subStorm condition. That is. it assumes that Vg 
is negative, so that electrons are repelled and ions attracted. In these equations 
Vg Is an algebraic quantity, that is, the Sign is implicit. 
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Figure R. 1-D Model Space Substorm Environment 



Afi Id evident from the sketch Irt Figure <i, secohdary electrons due to ion 
Impact are not accounted for in this model. This Is because the intent here is to 
use material ch&rglng characteristics found fitting the groi nd test model to 
®5tperlmental data In conjunction with the space model to predict space (charging. 
Since there were no ions used in the experiments, no coefficient for secondary 
electrons due to ion impact was determined. Therefore this current density 
Source is not considered. 'The surface voltage values predicted by this model are 
therefore somewhat largtr than if secondary electrons due to ion impact had been 
included. For example, a Secondary coefficient of 1 would yield about a 10 percent 
reduction in the equilibrium voltage calculated for aluminum at - r> kV. 

The procedure used to calculate charging is identical to that described for the 
ground test model. Material characteristics used were those determined by fitting 
the ground test model to the data. The capacitance of the aluminum (considered to 
represent the spacecraft "ground") was taV to be 15 V 10*’^ farads, which is 
the capacitance of a one meter diameter sphere (to represent a "typical" space* 
craft dimension) to infinity. This capacitance was chosen because the relevant 
capacitance f. r charging floating metal objects is that of the object to its surround- 
ings. Those parameters (notably C and ? for Teflon, and K for aluminum) which 
Varied as functions of beam (and therefore surface) voltage were associated with 
the equilibrium surface voltage for the appropriate test for purposes of making the 
space Voltage calculations. The relationships botweeh electron and ion tempera- 
tures and between temperatures arid current densities were taken from the Pro- 
visional specification for the Geomagnetic Substorm Environment. This 
Specification is given as Figure 1. Thus, results of the space calculations, shown 
in Figure 8, are given aS functions of electron temperature only 

In Figure C(l?) two cuti'es are shown for the surface voltage of silvered Teflon 
as a function of electron temperature. The first curve calculated used the experi- 
mentally determined value of 9 V 10^'’ SNcm for the effective re.-^istivity of Teflon, 
■ftiis curve bends sharply to the right as electron temperature increases. An 
inspection of the current densities driving the equilibrium indicated that leakage 
current played a large part in the equilibration of the teflon, thus, this curve 
yields a "good" value for the surface voltage of Teflon if the spacecraft "ground" 
is actually near plasma ground. This would be the case if, for example, photo- 
emisFion were holding the spacecraft ground near plasma ground arid the Teflon 
surface of concern Vrere shaded. However, if the spacecraft is assumed to be in 
eclipse, the aluminum representing the spacecraft ground is predicted to acquire a 
large negative voltage (solid line). In this case it is dear that leakage < urrent 
can not drive the equilibration of the Teflon surfac«- according to the dashed line. 
Therefore a second curve Was calculated for the Teflon, based on the assumption 
that there was no leakage. This is shown In thO dash-dotted line. 
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f^igut^e 8(^) sh6ws th6 of Charging predicted for the ypace condition. The 
time to oharge the alufninuhi ts predicted to be eeoorida, comparable to the time 
required to charge the floating plate in the ground experimenta. The time required 
to charge the Teflon is predicted to be several minutes at lower electron tempera* 
tures, ranging to tens of minutes as the electron temperature increases. Thus, 
equilibration of Teflon is predicted to require Significantly longer time in space 
than it does in ground^experiments. 

Using the curves shown in Figure 8, it is possible to predict the response of 
a "spacecraft", composed of an aluminum structure partially covered with silvered 
Teflon,, to substorm and eclipse conditions. Such a set of predictions is shown in 
Figure 9. For purposes of this figure, it is assumed that photoemissioh Is Sufft* 
dent to hold illuminated surfaces close to plasma ground. 

Figure 9(a) assumes that the Teflon surface of Interest is Shaded when the 
Spacecraft iS in Sunlight. A Substorm injection With a 5 keV Maxwellian electron.... .. 
distribution is assumed to occur at t = 0 and this environment is assumed to 
remain constant throughout the time shown. An..eclipse is assumed to occur from 
t = 60 min to t = 120 min. When the SubStorm occurs, the dark Teflon charges 
according to the dashed curve of Figure 8(a), Shown in Figure 9 as a solid curve; 
the aluminum at "ground" is assumed to be held near plasma ground by photo- 
emission. ...The Teflon..Snrfari^..-re.arhjes..itg leakage dominated equilibrium voltage 



Figure 9. Response to Substorm and Eclipse. 1-D Model 
Predictions 
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of -fl Ti kV with Itf time constant of about 20 min. When the epacet raft cntei’e 
eclipse, the aluminum charges oulckly (in seconds) to Its equilibrium voltage of 
- 10 kV It is assumed that the charge on the Teflon surface Is immobile oil this 
time scale, so this surface remains at -8. 5 kV. At this point, the leakage current, 
which was driven by the voltage differential of 8. 5 kV disappears, so the Teflon 
finds Itself to be no longer in equilibrium with Its environment and proceeds to 
charge to its "floating" value of —14. 5 kV in a characteristic 20-30 mlrt period. 
Upon exit from eclipse, a similar pattern is followed, the aluminum falls quickly 
to near plasma ground. Because this discharging IS driven by phOtoemlssion. it 
requires only about 0. 02 sec for the aluminum to roach plasma ground (assuming 
-10‘® A/cm^ photo current). Now. the Teflon again finds Itself out of equilibrium 
with its plasma environment, and-proceeds to discharge slowly to its previous 


equilibrium potential of -8. 5 kV, 

Figure 9(b) Shows a similar type of time history for an Insulating surface 
which IS exposed to sunlight. Again, the solid line represents the surface voltage 
of the Teflon and the dashed line the spacecraft ground. The entry into eclipse 
and subsequent charging up is analogous to the charging of the composite sample 
discussed in Section 4. The aluminum charges rapidly (in seconds) to its equi- 
librium value. Because the Teflon had ho significant charge on its surface, its 
voltage follows that of the aluminum until the aluminum reaches equilibrium. The 
Teflon then continues to charge slowly to US equilibrium potential. Upon exit from 
eclipse, both the aluminum and the Teflon are discharged by photoemisSion. Fhus, 
the aluminum reaches plasma ground in about 0. 02 sec, as in the previous rase. 
The Teflon also discharges more quickly than it charged: it requires about 4 min 


to reach plasma groutid. 

These results indicate a need for charging studies which take into account 
relative charging rates as well as different equilibrium charging levels of various 
spacecraft surfaces. A "typical" spacecraft has several different types of sur- 
faces (solar cells, thermal blankets, etc. ) each of which can be expected to charge 
with its own time constant. The importance of the effect of the different time 
constants should be assessed. 


. 1 . (.oM.i.i irtMi iu.'i vrth> 

The present study has resulted in the development of a set of "material charg- 
ing characteristics" which describe the charging of small (300 cm ) samples of 
5 mil silvered Teflon and oxidised aluminum. Based on these characteristics, 
predictions of charging in space have been made and used to estimate the behavior 
of a composite body under conditions of s ibstorm and eclipse. Several interesting 
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dlffrrencre charging behavior Under ^?round teat rondltions and predii ted 

behavior in 8poc^ have been noted foe silvered Terton. Under ground te^t c ondi* 
tions Bilveeed T*efton acquire the 8ame surface potential fhourited Oh «roun‘a d and 
floating substrates because the equilibration is domltirted by surfnee emission 
phenomena. However, in space, significant differences arc predicted in surface 
voltage for those two ihountirig configurations, tliis is because leakage current 
dominates the equilibration of the Teflon When the aluminum is at ground, while 
surface emission dominates for the floating aluminum case. The time required 
for silvered Teflon to charge to equilibrium irt tests is several minutes; this time 
is predicted to be several tens of minutes ih space. 

In contract to ttie several minutes to several tens of minutes time scales for 
silvered Teflon, floe ting aluminum samples arc observed in ground test and pre- 
dicted in space to charge to equilibrium in seconds Their capacitance is muc h 
lower than that of Teflon since it is determined by their surroundings. There are 
also orders of magnitude differences in time scales for discharging by photo - 
emission of aluminum and Teflon. This discrepancy in charging and discharging 
rates gives rise to sudden changes in the electric fields which the Teflon must 
sustain upon entry into and exit from eclipse. It is felt that these differential 
charging rates as well as differential charging levels may be important and should 
be investigated further. Thus, transient (quesistatic) as well as steady state 
models should ae developed for charging. 

The One-dimensional models described herein have been found useful in the 
interpretation of experimental results, and as guides to relating tesi results to 
expected space behavior. Models of ground test situations are needed since they 
can be used interactively with test data. This is especially true since it is im- 
possible, or at least impractical, to simulate accurately the geosynchronous 
environment. Thus environment scaling must be done through use of models, at 
least for the present. 

t'irtally, higher dimensional models are needed. One- dimensional models 
can not account for such things as edge effects or interactions between adjacent 
surfaces at different potentials as with different charging properties. Such effects 
are clearly important,^ and may dominate the charging behavior of multisurface 
samples and spacecraft, 
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Api^ndU A 

Ohe’dimenaiOAal Ground Test Model 


I. m4tUWXWUS 

In this modf’l,-. t'lef'trons from tlu* iAoviton gun nro assumed to approach the* 
planar sample normally. All motion Is restricted to the x direction (s(‘(‘ Figure 
Al)» The electron beam is assumed monoenergetic, with energy 


where e Is the electronic charge and V the beam voltage. The current density 
emitted from the gun is given by 



where n^ = particle density. It is assumed Viere that n^e (the charge density) is 
constant, in order to account for the spreading of ttie beam in the real situation. 
Thus the continuity equation requires that some particles are ’Tost”. 

We wish to calculate current densities to the sample surface. Current densi- 
ties to be considered are those due to primary electrons, secondary electrons, 
backscattered electrons, and leakage through the bulk Of the insulator. The insula* 
tor is assumed to be mounted above a grounded substrate for purposes of calcula- 
ting leakage. 

Throughout this dCvel ^pment the sample surface is assumed negative, and all 
secondary and backscattered electrons are assumed to escape. All signs are given 
explicitly so that symbols represent positive quantities. 


f^ROY 
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Figure A 1. 
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I I t hhi:M iilA^m iH i: TO iMmuin n,i:i 


C’onpldrr eleotronfi appropohlng tht* Ramplc- with onfr^y i: It Ifi riSButmcl that 

the cleotronp do not roUldt' with onf anoth< r, and that motion is la otriftrd to ono- 
dimenslon. If th«* sample surfoee has a repulsive potential of maf;nitude Vg, e nergy 
eonservation requires that they arrive at the surfaer- with energy Sim e 

all eleetrons leave the gun with energy r>V|j. the eurrent density to the sample sur- 
faee is simply 


/ 2 \1 1 '' 

J<. • "o- iTS-) 


Rt^arranging and using Kq. (Al), wr* fiiHi 




since we are requiring that n^e remain constant. 


3. C( ItIttAt IIKNSih III K TO Sl:i:ONO\I(\ CI.KCTIUINS 


Sternglass* has given the following expression for secondary yield as a func- 
tion of primary electron energy at impact: 


6(E.) = 7.46 

i m eV. 




where is the maximum yield, is the primary energy for which maximum 


yield is attained and E. is primary energy at impact From the discussion given 


in Eq, (A2) above, for this case 


El = . .Vg 


The secondary cur 'ent density is then 








r«i- •' 


which is 


' 

i' 

■VI 

% 

1: 



(A7) 


4. CIJHKKNT DENSITY DUE TO BACKSCATTEREl) ELECTRONS 

No analytical oxpresalon was found for backScattered electron emission. Fat 
simplicity, it vtfaS therefore assumed that backscattered electron current density 
represents a fraction of the Incident current density. Thus, 

/ 





5. LEAKAGE CURRENT DENSITY 

Leakage Current is generally represented by 

, .Zs. (AS) 

n ■ R 

In terms of bulk resistivity this iS 

, .1^ (AlO 

h~ fit - 

Miihere .? is.bulk.,reSiS.tlV^y, A IS the area and the thickness of the Sample. Then 

1 -JL- Zs (All 

M ■ A " 

6. THE 1-D MODEL 

The piimary electron current density i‘epresents a source of electrons arrlv 
Irtg at the sample, the other three current densities represent loss of electrons 
from the surface, thus the net current density to the surface IS 


M 
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ic ' " ^BS * 




This net current density plays the role Of a Changing curi'ent to the surface, thus, 
if we represent the sample's charging as the charging of a capacitor, we have 


dV,j 

Jc ' ^ dt “ ’ ^BS ' ^e 


(A 13) 


where C is capacitance, here expressed in farads per stjuare centimeter to main*- 
tain consistency of units. Equation (A 13) iS solved in the manner described in the- 
test on a compUter-to calculate the charging. 

It remains to associate the experimentally measured parameters. Surface 
voltage, and total current to ground, with calculated values. The surface voltage 
association Is trivial; it is Simply -Vg. The total current to ground is the charging 
current plus the leakage current for the sample as a whole. Thus 




A(3c + ii> 


(A 14) 
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Appendix B 

()n6-Dim0nSi«nol Spdce SubSMrm Modil 


1. INTHOIHICTION 


The one-dimensional space model assumes a two-dimensional isotropic 
Maxwellian Velocity distribution tor primary particles, and a spherical collection 
geometry. The calculations are essentially those for a spherical Langmuir probe. 
The present calculations are based on the work. of Langmuir* and Grard.et al. ^ 
and follow closely the derivation of Cauffman. ^ The latter work has not been 
published; therefore portions of it ar e reproduced here for clarity. Such portions 
are identified by superscript reference. 

Geometry for the calculations iS depicted in Figure Bl. ^ The sheath is as- 
sumed to have radius a, and the collector radius R. The radial and tangential 
velocity components in the "undisturbed" region (sheath edge) are v^. and v*. res- 
pectively. and those at the collector Surface u^ and u^. The Surtace potential of 
the collector is Vg. The potential In the ShSath is assumed to be a function of 
radial distance from the collector and to be monotOnic. The plasma is assumed 
cOlhsionlesS, that is, orbit limited theory applies, and energy and angular mo- 
mentum are assumed constant for each particle. 

"nie integral requiring solution for current densities due to primary elec- 
trons •' and ions, and backscattered electrons is 



Figure Bl. 
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(m) 


J. = 


4irR^ 


4ff a 
4»R^ 


J. 




R dE 


R 


( Rr - O ,,) 


wh^^6 the subscript I startds for bither o (elfectrortS) or p (H+ ions), and 0^^ and Ej^ 
are the angle and energy at impact on the coUectoJ* (see Figure Bl). 

For Secondary electron current, solution must be found for ' 



2 M 

f f 


(B2) 


where s (E„) is the secondary electron yield aS a function of electron impact 
K 

energy. 

These integrations cannot be performed dir.ectly because the distribution 
function for. the particles at the collector is unknown, and therefore we can not 
determine dj^/dEp. However, we do know the distribution at the Sheath edge, and 
can therefore determine dj./dE^. If we assume the plaSiiia to be colllSiOnleSS, we 
can also convert the limits on Ep and 0p to limits on E^ and and perform the 
required integrations On theSe variables. 

l.l i'ohdilions for t'ollwlion 

In order to contribute to current collected at R, a particle must have energy 
> 0 and direction Op £ 6 s ir/2. Since the plaSma is assumed colliSlonless 
and V(r) is assumed monOtonlc, each particle's energy and angular momentum 
must be conserved. Assume the particles of interest have charge -e. Energy 
c .nServatlon demands 


lmj(u| + = I m. fVj? + v^^) + eVj^ Ep 


" E^ + eV, 


Angular momentum conservation demands 


(B3 


RUj = avj RU sin 0p = av Sin 0^ 


Mtfhere 


r 2., ,2.1/2 

u = (u„ ) 


and 


/ 2 . 2 , 

= (v + Vj ) 


1/2 


(B4 


/ 


Solving fol* and 6^ Ih tbrms of Ep and fljj yields the conditions for collection^ 




» > Ep £ 0 


-1 A 


« =* E sE 
a o 


^ 0 for Vg > 0 (attraction) 

-eVg for Vg £ 0 (repulsion) 


1.2 BOcrity Format 

An isotropic Maxwellian v.elocity distribution in three -dimensions has a 
distribution given by 


2^J 


-m^v ' V 


We are interested ih a two-dimensional distribution which can be found from 


g(v^. v^) = / f(v^, ’ ^ Sin \ , V^ cos X ) d\ 


Substituting for f(v) and integrating, we have 


, 1 /m. V*'" |-m^2 + v2) 

Langmuir^ gives the incremental current across the Sheath as 


d. = 4 »f a^ n e v^ g(Vj., Vj) dv^. dv^ . 


Then. 


.i Ih + v,2)l 

^ [- Te f - 


arid, chaflgirig to E , § coordinates We Kav^"^ 

a & 


i 

f 


.s’* 


<eV^r 


E^ exp 


[•^J 


sin cos d0„ dE 
a a a s 


(BID 


as required. 


2. I^HIMMIV AM) ll\i:KSC\TTtlU:i) PAKTICLK Cl RKUNT DKNSITICS 


Since backScattered electrori current density is considered to be simply a 
fraction of the incident electron current density, the Same integration applies to 
both. The calculation for positive ion collection is the Same as for electron Collec- 
tion with appropriate sign changes to account for the positive charge, and using the 
ion mass and temperature. In Eqs. (B9) to (Bll), the sign on the leading charge 
has been Suppressed. Appropriate sighs will be supplied in Section B. 5. 

The integration to be performed is 




(. V ,)2 i J 


sin 9 cos 6 
a c 


exp 


(■*;) 


Ids dE . 
a a 


(B12) 


3 1 

Cauffman evaluates this integral and finds, in agreement with Langmuir 



Vg s 0 (Repulsive) 
Vg > 0 (Attractive) 


(B13) 


where 

/eVi j/2 

■ ”i® \2»m^ 

Since the interest here is in modeling charging In substorms with nophotoemiseion, 
we ejcp^ct Vg negative. Thus, electrons are repelled ahd ions attracted. So we 
have for electrons 


3. 


eo 
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(S14) 




tot iotid 


•J- 




ip-t 


'po(-- 




<B15) 


Where the minus Sign reflects ths fact that ihht' are attracted by negative Vg. For 
backscattered electrons. 


Jbs = *^3eo 


(Blfi) 


where S is the backscatter coefficient. 


3. SECONDARY ELECTRONS DUE TO ELECTRON IMPACT 


Secondary yield aS.a-function of electron impact energy has been given by 
SterngiabS^ aS 


6{Ep) = 7.4 6 


m eVZ ®*P 
m 


r\ 


, 1/2 




(A4) 


where 6^ is the maximum yield and eV^ the energy at which the maximum yield 


is Obtained. To determine secondary electron current density, we must multiply 
the left hand side of equation (Bii) by 6 (Ej^) and integrate. Thus we need 




2 m 


S ° (eV^ ®''m j. 





X exp 




/2 


■iV^ ‘2 V"^. 


sin cos dfl^ dE 
e a a a 


(B17) 


Which is Cauffman's^ equation for secondary electrons, except that he uses a 
sec dependence of t on 9^, which is not used here. The 0^ integral is the Same 
as before and yields 


j B j 


7.46^ r ^ 


(eV„)^eV^ 


-®a J 



J 


dE. 


(018) 
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) 




and subfld- 


Now. seqU6ntlally setting + eVg 

tuti'ig, we find 


and n = (eVg)'^^^ 





ejip i-n^) dn 


fBig) 


where 

'’o 



Since we are considering substerm eases only here, the condition Vg £ 0 Is of 
Interest. For this case the Integral In Eq. (B19) IS just 



(B20) 


Recalling the expression for repeated Integrals 


of th^ error function complement 


5 


erfc(fc) = 



<t-z)" 

n* 


6xp (“t^) dt 


(B21) 


where, by definition. 

o6 

l" erfe(z) = erfc(t)dt 

t 

and 



exp (-t^)dt “ ei*fc(a) 




and Identifying, t, z and h with the afftjrofirlate varlablee In (B20), we have 



So, from Etjs, (B19) and (022), we have 



and we note that the dependence of jg on Vg is the same as that of and j^g. 


( 022 ) 


(02 3) 


1. .LKAkAliK CntKKNt DKNSITV 

Leakage current density is defined in the same manner for the space model as 
it was for the test facility model (see Se-tion 5 of Appendix A). Thus, we have 



where p is bulk resistivity and 1 is the thickness of the insulating film. 

3. TIIK Ml SP \I;K MOtlW. 

Now, the net current density to the sample surface is 

4 = + i +1 + i„_ - i, (025) 

^c ^p •’S •'BS '^1 

where the Signs on the current densities are given explicity here, and we recall 
that Vg in this model Is algebraic (that is, can be positive or negative), although 
the derivations have assumed it negative. 

The net current density plays the role of a charging current to the surface, so 

that 


i = C 
•’c 


'‘"i- i 

‘dt ■ 


e ^s ^BS ” ^1 


(02 f!) 
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3^ 


mm 


wheft C is capacity expressed ift fal*adfi per square centlrhctor and where we 
assume we are chargi .g a capacitor. Equation (B2fi) is solved in the same way 
as Eq, (A 13) on a computer to determine Vg versus time for charging. 
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